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Microgrid  is  an  effective  concept  for  distributed  renewable  energies  connecting  to  utility  grid.  However, 
enormous  challenges  have  also  been  brought  into  the  microgrids  study  considering  the  intermittence, 
randomness  characteristics  of  renewable  energies.  As  a  result,  the  power  generations  of  distributed  gener¬ 
ators  fluctuate  frequently.  Furthermore,  it  is  difficult  to  obtain  the  parameters  of  microgrids  accurately. 
Thus,  the  study  on  small  signal  stability  of  microgrids  is  especially  important.  This  paper  presents  a 
review  of  exiting  small  signal  stability  methods  for  microgrids.  Besides,  a  new  approach  based  on  Sin¬ 
gular  Entropy  and  Matrix  Pencil  is  proposed.  Meanwhile,  its  validity  and  effectiveness  are  confirmed  by 
some  examples  compared  with  Prony  algorithm.  Then,  the  proposed  method  is  also  applied  into  a  micro¬ 
grid  with  an  improved  droop  control  strategy.  Finally,  a  framework  based  on  the  proposed  method  for 
small  signal  stability,  warning  and  emergency  control  online  in  microgrids  is  available. 
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1.  Introduction 

Recently,  the  shortage  of  fossil  fuels,  vulnerability  of  large  grid 
and  the  power  blackouts  caused  by  it  have  been  more  and  more 
serious  [1-3].  Therefore,  distributed  and  renewable  energies  such 
as  wind,  solar  and  so  on  gained  more  and  more  interests  [4-6]. 
However,  because  of  their  intermittence,  randomness  and  uncer- 


*  Corresponding  author.  Tel.:  +86  13588846066. 

E-mail  addresses:  zengerzheng@zju.edu.cn  (Z.  Zeng),  yanghuan@zju.edu.cn 
(H.  Yang),  rongxiang@zju.edu.cn  (R.  Zhao). 

1  Tel.: +86  13567123512. 

2  Tel.: +86  13906508946. 

1364-0321  /$  -  see  front  matter  ©  201 1  Elsevier  Ltd.  All  rights  reserved, 
doi:  1 0.1 01 6/j.rser.201 1 .07.069 


tainty  caused  by  meteorological  factors,  it  is  hard  to  connect 
renewable  energy  sources  to  utility  grid  directly  [7-10].  By  inte¬ 
grating  distributed  and  renewable  sources,  energy  storage  devices, 
a  variety  of  loads,  data  acquisition  and  supervisory  control  devices, 
microgrids  realize  the  interface  between  the  distributed  renewable 
sources  and  utility  grid  [11-13]. 

At  present,  the  researches  on  microgrids  mainly  focus  on  the 
optimal  design  of  system  capacity  [14],  power  generation  opera¬ 
tion  and  control  [15-18],  system  modeling  and  simulation  [19-23], 
power  quality  [24,25]  and  so  on,  while  rarely  on  the  stability 
[26,27].  However,  stable  and  reliable  power  supply  to  loads  or  util¬ 
ity  grid  is  the  basic  function  of  microgrids.  So,  study  on  stability 
of  microgrids  is  very  necessary  and  important.  Unlike  transmission 
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Nomenclature 

P;.  Qi 

active  power  or  reactive  power  output  of  the 
ith  distributed  generator  (DG),  DG/,  whose  output 
apparent  power  is  S/  =  P/  +jQ/.  Especially,  Pre^  (Qre/z) 
is  its  reference  value.  Similarly,  Pg  ( Qg )  and  Sg  are  the 
ones  supplied  by  utility  grid 

P,  Q 

active  power  or  reactive  power  after  digital  filter  in 
inverter  controller 

hi>  Qli 

active  or  reactive  power  demand  of  the  ith  load 

kpp*  kpi 

proportion  and  integration  coefficient  of  DGZ  for 
active  power  droop  control 

kqp*  kqi 

proportion  and  integration  coefficient  of  DGZ  for 
reactive  power  droop  control 

r,  xL 

resistance  and  inductance  of  line  per  kilometer 

co  o 

rated  angular  frequency  of  the  microgrid 

Eo 

rated  potential  of  DG 

Kef 

three-phase  reference  voltage  vector  of  DG 

Ep 

potential  of  DG  injected  into  the  microgrid 

Ts 

sample  time  of  raw  data  or  sampled  data  for  small 
signal  stability  modes  extraction 

y 

raw  data  acquisited  from  microgrid,  and  y(i)  (i  =  1 , 2, 

. . .,  N)  is  its  ith  point,  where  N  is  the  length  of  it 

Y 

formed  Hankel  matrix  by  raw  data 

Ek,  A Ekt 

A E'k  the  kth  Singular  Entropy,  Singular  Entropy 
Increment  and  Differential  of  Singular  Entropy 
Increment,  respectively 

n 

modal  order  of  raw  data 

Hi 

the  eigenvalue,  frequency  and  damping  ratio  of  ith 
mode 

SNR 

signal  noise  ratio 

network  or  distribution  network,  microgrids  can  operate  in  isolated 
mode  or  interconnected  mode  [  1 2 ].  This  feature  brings  out  a  serious 
challenge  to  stability  analysis  of  microgrids.  Stability  of  microgrids 
also  includes  two  parts.  One  is  small  signal  stability  (SSS)  in  small 
disturbances  sense,  and  the  another  is  the  transient  stability  in  large 
disturbances  sense.  Except  for  large  disturbances,  the  power  out¬ 
puts  of  wind  turbines,  photovoltaic  arrays  and  other  renewable 
sources  are  time-variant,  which  can  be  regarded  as  small  distur¬ 
bances.  What  is  more,  they  may  couple  with  load  fluctuations  too. 
Thus,  these  kinds  of  small  power  fluctuations  are  so  frequent  that 
they  affect  the  operation  status  of  microgrids  constantly.  In  a  word, 
the  difficulty  and  necessity  of  study  on  SSS  of  microgrids  is  evident. 

This  paper  reviews  the  current  approaches  for  SSS  of  micro¬ 
grid.  Besides,  an  improved  Matrix  Pencil  combined  with  Singular 
Entropy  is  presented  and  applied  to  extract  the  SSS  modes  of  micro¬ 
grids.  The  paper  is  organized  as  follows.  In  Section  2,  a  review  of 
exciting  and  prosiable  methods  for  SSS  of  microgrids  is  available. 
In  Section  3,  the  basic  mathematical  model  of  Singular  Entropy 
Matrix  Pencil  is  proposed.  In  addition,  the  microgrid  structure  and 
an  improved  control  strategy  for  DGs  are  provided.  As  well,  some 
simple  examples  are  done  to  confirm  the  proposed  method.  Then, 
the  proposed  method  is  applied  to  the  microgrid  for  its  SSS  modal 
parameters  extraction.  Section  4  is  the  discussion,  a  framework 
based  on  the  proposed  method  for  the  SSS  identification,  warn¬ 
ing  and  emergency  control  online  in  microgrids  is  presented.  The 
conclusions  are  summarized  in  Section  5. 

2.  Study  on  mall  signal  stability  of  microgrids 

Although  there  are  a  few  approaches  for  SSS  of  microgrids  so 
far,  their  disadvantages  are  obviously.  Therefore,  a  further  work 
is  strongly  demand  and  a  new  field  needs  to  be  explored.  In  [26] 


a  method  for  SSS  of  microgrids  in  interconnected  mode  had  been 
proposed.  A  model  presented  in  [27]  is  suitable  for  SSS  of  micro¬ 
grids  in  isolated  mode.  They  both  featured  the  eigenvalue  analysis 
approach  as  their  theoretical  basis  and  established  SSS  models  in 
interconnected  or  isolated  mode,  respectively.  However,  results  of 
eigenvalue  approach  for  a  fixed  system  may  deviate  far  from  the 
actual  condition,  due  to  the  perturbation  of  system  parameters, 
load  fluctuations,  power  scheduling,  distribution  network  distur¬ 
bances  and  other  causes  of  small  disturbances  in  the  real  world.  In 
addition,  the  calculation  process  of  eigenvalues  depends  on  system 
operation  status  strictly,  unfortunately,  which  is  hard  to  meet.  The 
system’s  operation  status  varies  frequently  and  largely,  since  the 
power  outputs  of  wind  turbines  and  photovoltaic  arrays  always 
vary  in  a  large  range  slowly.  As  a  result,  it  brings  tremendous 
challenges  to  off-line  SSS  analysis  based  on  eigenvalue  analysis 
approach.  What  is  more,  conventional  matrix  eigenvalue  approach 
is  based  on  state  space  matrix  obtained  by  the  local  linearization 
model  of  the  complex  nonlinear  system.  If  we  consider  sophisti¬ 
cated  control  strategies  as  well  as  a  number  of  nonlinear  couplings, 
the  eigenvalue  approach  will  not  be  a  good  solution.  Sometimes,  the 
linearization  model  even  may  not  be  available.  Besides,  the  inter¬ 
action  among  different  modes  may  cause  new  oscillation  modes 
based  on  normal  form  theory  [28-30].  Modes  interaction  is  abso¬ 
lutely  ignored  by  eigenvalue  approach,  because  it  is  just  the  results 
of  first-order  Taylor  expansion.  Especially,  its  results  of  SSS  arm  to 
the  whole  system.  As  a  result,  we  cannot  know  the  specific  modes 
for  each  DG  or  take  emergency  control  strategies  corresponding  to 
the  right  DGs,  when  there  are  unstable  modes  in  microgrids. 

As  mentioned  previously,  SSS  of  microgrids  may  be  calculated 
by  eigenvalue  theory  analytically,  but  there  are  some  serious  dis¬ 
advantages.  On  the  other  hand,  the  SSS  modal  parameters  of 
microgrids  can  also  be  identified  or  extracted  by  online  data 
directly,  which  can  overcome  the  disadvantage  of  eigenvalue  anal¬ 
ysis  approach.  This  kind  of  methods  has  a  wide  range  of  applications 
in  transmission  network  [31  ],  although  they  also  have  some  draw¬ 
backs.  Common  identification  methods  are  Prony  algorithm  [32], 
Hilbert-Huang  Transform  (HHT)  method  [33]  and  the  Matrix  Pen¬ 
cil  method  [34,35],  and  so  on.  Traditional  Prony  algorithm  uses 
complex  exponential  functions  to  fit  the  raw  data,  so  it  may  pro¬ 
duce  a  large  number  of  false  modes  beside  its  sensitivity  to  noise. 
Although  some  improved  iterative  Prony  algorithms  can  partially 
inhibit  the  effect  of  noise,  they  also  increase  the  computation 
complexity  [36-39].  HHT  based  modal  identification  method  can 
extract  the  instantaneous  modal  information,  but  it  results  to  a 
high  computational  complexity.  Furthermore,  its  end  effect  needs 
to  be  overcome  [40,41].  Matrix  Pencil  method  is  a  typical  signal 
recognition  method,  too.  Recently,  it  has  been  applied  to  power  sys¬ 
tem  low  frequency  oscillation  in  [42],  as  it  can  overcome  the  noise 
interference  and  have  low  computational  complexity  [34,35].  Tra¬ 
ditional  Matrix  Pencil  method  requires  a  threshold  to  obtain  the 
modal  order,  which  is  depends  on  whether  the  ratio  between  each 
singular  value  and  the  largest  one  bigger  than  the  threshold  or  not. 
How  to  set  the  threshold  is  affected  by  the  size  of  raw  data  and 
the  strength  of  noise,  which  is  hard  to  be  generalized  and  has  no 
appropriate  theoretical  guidance. 

A  lot  of  relevant  analysis  showed  that  the  key  of  the  modal 
parameters  extraction  from  raw  data  for  SSS  study  is  how  to  deter¬ 
mine  modal  order  accurately.  If  the  modal  order  is  few  than  its  real 
value,  it  will  miss  some  modes.  In  contrast,  if  the  modal  order  is 
more,  it  will  cause  over-fit,  which  results  to  false  modes  and  reduces 
the  accuracy  of  the  modal  extraction  results. 

In  this  paper,  Singular  Entropy  is  proposed  to  determine  the 
modal  order  of  the  raw  data.  It  can  not  only  obtain  the  modal  order 
whether  the  noise  affects  the  raw  data  or  not,  but  also  can  remove 
the  false  modes  and  reduce  computational  complexity  automati¬ 
cally.  This  method  also  uses  Matrix  Pencil  method  to  extract  the 
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modal  parameters  and  analyze  the  SSS  of  the  microgrids  based  on 
actual  operation  data  in  real  world.  So  the  SSS  information  of  micro¬ 
grids  can  be  obtained  real-time  and  online  dynamically.  It  can  adapt 
the  operation  status  and  small  signal  modes  varying  with  time  in 
microgrids.  On  the  other  hand,  it  can  meet  the  demand  for  security 
and  stability  warning  and  emergency  control. 


Increment,  it  is  easy  to  extract  the  order  of  the  corresponded  inflec¬ 
tion  point,  namely,  the  modal  order  n. 

3.1.3.  Form  a  new  singular  matrix 

After  we  determined  the  modal  order  n,  a  new  matrix  D'  can  be 
formed  by  the  first  nth  singular  values  of  D: 


3.  Singular  Entropy  Matrix  Pencil  method  and  its  microgrid 
application 

3.1.  Mathematical  model  of  Singular  Entropy  Matrix  Pencil 
method 

Firstly,  the  Singular  Entropy  is  presented  to  determine  the  order 
of  the  signal  modes,  which  can  accurately  determine  the  modal 
order  regardless  of  the  noise.  Then,  Matrix  Pencil  method  is  used  to 
extract  the  modal  parameters  of  raw  data.  Combined  with  Singular 
Entropy  and  Matrix  Pencil,  the  proposed  method  can  remove  the 
false  modes  and  reduce  computational  complexity,  automatically. 
The  process  can  be  described  as  follows. 

3.1.1.  Form  Hankel  matrix  structure 

The  acquired  or  simulated  raw  data y(i)  (i=  1,2,...,  N)  can  form 
a  Hankel  matrix: 

y(i)  y(  2)  y(i  +  i)' 

y(  2)  y(  3)  y(i  +  2) 

y=  :  :  :  (1) 

_y{N-L)  y(N  - L  +  1)  •••  y(N) 

where  L  is  the  matrix  pencil  parameter  and  a  proper  chosen  L 
can  suppress  noise  iny.  Take  L  =  N/4~N/3,  usually.  For  Y ,  we  have 
Y  =  UDVt  according  to  the  singular  value  decomposed  (SVD),  where 
singular  value  matrix  D  is  diagonal  Dg#^x^l+11  Its  element  a a 
(i  =  1 ,  2,  . . .,  m)  is  the  ith  singular  value  of  Y ,  here  m  =  min  {N  -  L, 
L  + 1 }  is  the  minimum  dimension  of  matrix  D.Ue  r(n~lMn~l)  is  the 
left  eigenvector  of  Y ,  while  V e  |?(L+1  )x(l+ 1 )  is  the  right  one. 


3.1.2.  Determine  the  modal  order 

Singular  Entropy  had  been  introduced  to  detect  fault  in  [43].  In 
this  paper,  we  used  it  to  extract  the  modal  order  of  raw  datay.  Now, 
define  the  Singular  Entropy  Increment: 


where  the  kth  order  Singular  Entropy  Ek  is  defined  as: 


k 

Ek  =  Ej  k  <  m  (3) 

i=i 

So  the  kth  order  Differential  of  Singular  Entropy  Increment  is: 


AC  _  AC 

AE'k=  ki(k_  1 ) 1  =  AEk  -  k  =  2,3,..  .,171  (4) 

where  A E'k  indicates  the  Differential  of  Singular  Entropy  Incre¬ 
ment  AEk.  When  the  effective  signal  saturates,  the  Singular  Entropy 
Increment  converges  to  a  bounded  value  rapidly.  So,  there  is  an 
inflection  point  which  is  corresponding  to  the  signal  modal  order. 
After  that  point,  the  Singular  Entropy  Increment  is  small,  which 
can  be  considered  as  the  effect  of  noise.  Due  to  this  advantage  of 
Singular  Entropy  Increment,  no  matter  how  much  the  noise  is,  the 
increment  of  Singular  Entropy  will  converge  to  a  bounded  value 
with  a  noticeable  jump  when  the  information  of  effective  signal 
tends  to  saturation.  With  the  help  of  Differential  of  Singular  Entropy 


G\  0  0  •  •  •  0 

0  o2  0  •  •  •  0 

D=  ...  :  (5) 

0  0  0  •  •  .  On 

_  ^(N-L-n)xn 

where  D'  eR^N^xn,  its  first  n  rows  are  mainly  formed  by  the  first 
nth  singular  values  of  D  and  the  last  N-L-n  rows  are  zeros.  The 
formed  matrix  D'  can  eliminate  noise  effectively,  because  it  has 
removed  the  singular  values  caused  by  noise. 


3.1.4.  Obtain  modal  damping  ratio  and  frequency 
Now,  form  two  matrixes  Y\  and  Y2: 

(6) 
(7) 


Y ,  =  UD’V \ 
Y2  =  UD’V  2 


where  Vj  is  formed  by  the  1  to  L  rows  of  the  right  eigenvector 
VT  corresponding  to  first  n  singular  values;  V2  is  formed  by  the 
2  to  L+l  rows  of  the  right  eigenvector  VT  corresponding  to  first 
n  singular  values.  Based  on  Y\  and  Y2,  a  matrix  pencil  Y2-XY\ 
can  be  derived.  Modal  parameters  are  included  of  the  eigenval¬ 
ues  of  matrix  G  =  Y^Y2,  where  Y|  is  the  pseudo-inverse  matrix 
for  Y\.  If  the  number  of  non-zero  eigenvalues  for  G  is  M,  namely, 
A,(i  =  1 ,  2,  . ..,  M;M<n),  the  corresponding  modal  damping  ratio 
and  frequency  of  each  mode  can  be  obtained  by  Eqs.  (8)  and  (9): 


_ — sgn{Re(Xi )) _ 

Vi  +(/m(In  A,-/rs)/Re(ln 
Re( In  a,/T5) 


(8) 

(9) 


where  Xj  is  the  ith  eigenvalue  of  ith  mode.  and  f  are  the  damping 
ratio  and  frequency  of  the  ith  mode,  respectively.  Ts  is  the  sample 
time.  According  to  Eq.  (8),  the  system  is  unstable  if  any  damping 
ratio  is  negative.  That  is  a  simple  judgement  for  SSS  of  microgrid. 


3.2.  Microgrid  structure  and  control  strategy 

The  structure  diagram  of  a  microgrid  for  study  in  this  paper  is 
shown  in  Fig.  1.  It  integrates  10  kW  photovoltaic  arrays  (DGt ),  3  kW 
micro-gas  turbine  (DG2)  and  2  kW  fuel  cells  (DG3).  There  are  some 
local  loads  at  each  DG’s  terminal,  which  are  Load  1 ,  Load  2  and  Load 
3,  respectively.  DGs  are  connected  to  the  point  of  common  coupling 
(PCC)  by  lines,  where  load  4  attached.  Microgrid  connected  to  utility 
grid  at  PCC  simultaneously.  The  utility  grid  is  a  low  voltage  distribu¬ 
tion  network,  whose  rated  root  mean  square  value  of  phase-ground 
voltage  is  220  V  and  its  rated  frequency  is  50  Hz. 

Conventional  droop  control  of  DG  is  shown  in  Fig.  2.  According 
to  the  active  and  reactive  power  droop  characteristics  of  DG,  its 
decoupling  control  feature  of  frequency  and  potential  amplitude 
can  be  achieved  [44,45].  As  shown  in  Fig.  2,  the  conventional  droop 
control  strategy  can  be  expressed  as: 

(  co  =  coq  —  kp(P  —  Pref )  m  n\ 

\E  =  E0-kv(Q-Qref)  1  j 

where  to  and  E  are  angular  frequency  and  potential  amplitude  of 
the  inverter,  respectively.  Meanwhile,  too  and  E0  are  the  ones  cor¬ 
respond  to  the  non-load  situation.  P  and  Qare  the  measured  output 
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Fig.  2.  Traditional  droop  control  strategy  of  DG.  (a)  Active  droop  characteristic 
for  frequency  control,  (b)  Reactive  droop  characteristic  for  potential  amplitude 
control. 


and  reactive  power  after  digital  filter,  respectively.  Vref  is  the  refer¬ 
ence  voltage  output  of  DG,  which  is  given  by 

'  Vrefa  =  E  sin  (cot) 

Vreffj  =  E  sm{cot  -  2tt/3)  (11) 

Vrefc  =  E  sin  {cot  +  2tt/3) 

In  Fig.  3,  Ep  is  the  potential  of  DG  injected  into  the  microgrid. 
Prefi  and  Qrefi  are  reference  power  outputs  of  DG*,  which  is  decided 
by  the  steady  operation  condition.  Gi(s)  and  G2(s)  are  the  transfer 
functions  of  Proportional  Integral  (PI)  regulators.  In  this  paper,  we 
take  the  follow  expression: 


active  and  reactive  power  components  of  the  inverter.  Pref  and  Qref 
are  the  reference  ones.  kp  and  kv  are  the  droop  coefficients  for  the 
active  and  the  reactive  power. 

The  conventional  droop  control  strategy  works  and  coupling 
control  can  be  achieved,  when  the  line  is  inductive  purely.  How¬ 
ever,  if  we  take  into  account  the  resistive  components,  the  angular 
frequency  of  DG  will  not  only  relative  to  the  active  power,  but  also 
reactive  power  [46].  These  resistive  components  may  be  line  resis¬ 
tance,  the  equivalent  resistance  of  the  filter  inductance,  parasitic 
resistance  of  the  power  electronic  components  and  so  on.  Sim¬ 
ilar  results  can  be  reached  for  the  potential  amplitude  [47].  On 
this  situation,  it  is  hard  to  decouple  the  complex  and  nonlinear 
system  by  the  conventional  droop  control  strategy  as  mentioned 
previously. 

From  (10),  it  is  can  be  observed  that  the  traditional  droop  con¬ 
trol  can  be  regarded  as  a  Proportional  Integral  Derivative  (PID)  type 
controller,  or  just  a  proportional  controller  exactly,  whose  input  is 
the  active/reactive  power  deviations.  In  order  to  achieve  the  decou¬ 
pling  control  of  the  microgrid  with  resistive  components,  we  add 
an  active,  reactive  power  control  cross  and  an  integral  component 
and  the  conventional  droop  controller  or  proportional  controller 
[46,47].  Then,  an  improved  droop  control  strategy  is  proposed,  as 
shown  in  Fig.  3,  where  cof  is  the  angular  frequency  for  the  filter 
whose  transfer  function  is  G(s)  =  ojfl{s  +  cof).  P  and  Q.  are  the  active 


kni 

Gi(s)  =  kpp  + 


G2(s)  —  kqp  + 


s 

kqi 


(12) 


According  to  Fig.  3(a)  and  Eq.  (12),  the  improved  droop  control 
strategy  can  be  written  as: 

CO  =  COq  +  kpp(Pfe/  “1“  Epi  f  ( Eref  P)dt 

Jo 

+  kqpiQref  ~  Q)  T  kqf  f  {Qref  —  Q)dt 

Jo 

E  =  Eq  +  kqp ( Qref  ~  QJ  +  kqf  f  {Qref  ~  QJjdt 

Jo 

Epp{Pref  ~  P)  +  kpi  f  {Eref  ~  P)dt 

Jo 


< 


(13) 


The  important  parameters  for  the  microgrid  and  controller  are 
shown  in  Table  1.  Notation  that  r  and  xL  are  the  parameters  of 
lines,  and  the  length  of  lines  are  shown  in  Fig.  1.  The  micro¬ 
grid  receives  Sg  =  2+j0.6kVA  power  from  utility  grid  in  the  given 
steady-state  condition  according  to  power  flow  analysis.  All  DGs 
take  the  improved  droop  control  strategy  proposed. 


Table  1 

Parameters  of  the  microgrid  system. 


Parameters 

Pref\ 

Qref 1 

Prep. 

Qre/2 

Prep 

QreP 

PS 

Qg 

Values 

10 1<W 

0.4  kvar 

3  kW 

0.4  kvar 

2  kW 

0.3  kvar 

2  kW 

0.6  kvar 

Parameters 

Pl  i 

Qii 

Pl2 

0l2 

PL3 

Ql3 

PL4 

Ql4 

Values 

5  kW 

0.1  kvar 

1  kW 

0.5  kvar 

1 1<W 

0.1  kvar 

lOkW 

1  kvar 

Parameters 

kpp 

kpi 

kqp 

kqi  r 

*L 

coo 

E0 

Values 

0.0005 

0.0001 

0.0001 

0  0.243  Q/km 

0.83  Q/km 

314rad/s 

220V 
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Fig.  3.  The  improved  control  strategy  and  structure  of  DG.  (a)  Block  diagram  for  control  strategy,  (b)  The  control  structure  for  hardware  application  of  DG. 


3.3.  Simulation  and  results 
3.3.1.  Some  simple  examples 

Now,  we  consider  some  simple  examples,  which  can  confirm 
the  validity  and  effectiveness  of  the  proposed  method.  If  the  active 


power  P(t)  in  per  unit  satisfies  the  following  relationship: 

P(t)=y(t)  =  y^t)  +y2{t)  +y3(0  +  e(t)  (14) 

where  s(t)  is  white  noise.  And  the  sample  time  is  Ts  =  0.01  s.  Three 
oscillation  signals  with  damping  are  expressed  as: 

y-l(t)  =  e~0A5t  sin(20  x  27 it)  (15) 

y2(t)  =  e~035t  sin(2  x  27 rt)  (16) 

y3(t)  =  0.2e_a25f  sin(0.5  x  27it)  (17) 


There  are  three  oscillation  signals  in  P(t),  and  each  signal  corre¬ 
sponds  to  the  two  modes  in  complex  frame.  So  the  accurate  modal 
order  is  2  x  3  =  6.  When  we  consider  the  signal  with  noise,  it  is  sup¬ 
posed  that  the  white  noise  sequence  meets  normal  distribution 
s(t)  ~  N( 0,  a).  Fig.  4  shows  the  Singular  Entropy  Increment  in  case  of 
cr  =  0,cr  =  0.025  and  o  =  0.05,  where  “o”  represents  the  turning  point 
corresponding  to  the  modal  order  of  raw  data.  From  Fig.  4,  it  is 
clear  that  the  Singular  Entropy  Increment  reflects  the  information 


of  modal  order  accurately  (both  are  6).  When  we  consider  the  sig¬ 
nal  without  noise,  the  Singular  Entropy  Increment  is  zero  after  the 
modal  order  n.  When  we  consider  noise,  the  Singular  Entropy  Incre¬ 
ment  also  has  an  apparent  transition  near  the  modal  order  n.  After 
n,  we  regard  these  A Ek  is  caused  by  noise,  which  is  very  small  and 
closed  to  a  bounded  value,  approximately.  So  we  can  use  Singular 
Entropy  Increment  to  identify  the  modal  order  of  raw  signal. 

3.3. 1.1.  Case  1:  raw  data  without  noise.  In  this  part,  we  consider  the 
case  that  raw  data  is  not  polluted  by  noise,  and  give  the  modal 
parameters  extraction  results  based  on  Singular  Entropy  Matrix 
Pencil.  Fig.  5(a)  gives  the  Singular  Entropy  Increment  and  its  dif¬ 
ferential. 

The  modal  order  result  from  Singular  Entropy  Matrix  Pencil  is  6, 
and  the  signal  noise  ratio  (SNR)  for  the  fitting  is  SNR  =  585.68  dB. 
In  the  same  case,  we  use  the  Prony  algorithm  with  6  orders 
as  a  comparison.  The  results  show  that  the  signal  noise  ratio  is 
SNR  =  151.19  dB.  Table  2  and  Fig.  6  show  the  results  in  detail. 

From  Table  2  and  Fig.  6,  it  is  easy  to  find  that  the  singular  Entropy 
Matrix  Pencil  method  can  identify  the  modes  of  original  signal  well. 
In  addition,  it  is  more  accurate  than  Prony  algorithm.  Especially, 
the  phase  result  of  Prony  algorithm  for  mode  1  is  far  away  from  the 
accurate  value,  namely,  the  bold  value  as  shown  in  Table  2. 


Fig.  4.  Relationship  between  Singular  Entropy  Increment  distribution  and  the  modal  order  of  raw  data,  (a)  Raw  data  without  noise  in  the  time  domain,  (b)  Differential  of 
Singular  Entropy  Increment  AF.  (c)  Singular  Entropy  Increment  A £. 
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Fig.  5.  Singular  Entropy  Increment  and  its  differential  to  obtain  the  modal  order,  (a)  The  results  for  raw  data  without  noise,  (b)  The  results  for  raw  data  with  noise. 


Table  2 

A  comparison  of  identification  results  for  raw  data  without  noise. 


Frequency  (Hz) 

Damping  ratio 

Amplitude  (p.u.) 

Phase  (rad) 

Real 

Identified 

Real 

Identified 

Real 

Identified 

Real 

Identified 

Prony  algorithm 

1  20 

20.0001 

0.0036 

0.00358 

1 

1.0045 

1.5708 

0.3142 

2  2 

2.00078 

0.0278 

0.0278 

1 

1.0035 

1.5708 

1.4451 

3  0.5 

0.5016 

0.0793 

0.0793 

0.2 

0.2005 

1.5708 

1.5394 

Singular  Entropy  Matrix  Pencil 

1  20 

20 

0.0036 

0.00358 

1 

1 

1.5708 

1.5708 

2  2 

2 

0.0278 

0.02784 

1 

1 

1.5708 

1.5708 

3  0.5 

0.5 

0.0793 

0.07933 

0.2 

0.2 

1.5708 

1.5708 

(c) 

1.5 

1 

.  0.5 

a 

^  0 


-0.5 


-1 


Fig.  6.  The  identification  results  for  raw  data  without  noise,  (a)  Comparison  between  raw  data  and  fitting  data  by  Prony  algorithm,  (b)  Error  for  Prony  algorithm  (R-F).  (c) 
Comparison  between  raw  data  and  fitting  data  by  Singular  Entropy  Matrix  Pencil  method,  (d)  Error  for  Singular  Entropy  Matrix  Pencil  method  (R-F).  R,  raw  data;  F,  fitting 
data. 
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Fig.  7.  The  identification  results  for  raw  data  with  noise,  (a)  Comparison  between  raw  data  and  fitting  data  by  Prony  algorithm,  (b)  Error  for  Prony  algorithm  (R-F).  (c) 
Comparison  between  raw  data  and  fitting  data  by  Singular  Entropy  Matrix  Pencil  method,  (d)  Error  for  Singular  Entropy  Matrix  Pencil  method  (R-F).  R,  raw  data;  F,  fitting 
data. 


When  the  raw  data  has  no  noise,  the  extracted  modal  parame¬ 
ters  can  meet  some  general  requirements  for  simulation,  such  as 
oscillation  modal  identification,  damping  ratio  analysis,  evaluation 
of  the  rationality  of  the  microgrid  system  structure,  control  equip¬ 
ment  capacity  and  location  planning  and  configuration,  optimal 
design  for  control  strategies  and  so  on. 

33 A 2.  Case  2:  raw  data  with  noise.  In  this  part,  we  consider  the 
situation  that  the  raw  data  in  (14)  has  noise.  We  give  modal  identi¬ 
fication  results  based  on  the  Singular  Entropy  Matrix  Pencil  method 
and  Prony  algorithm.  It  is  assumed  that  the  noise  meet  e(t)  ~N(0, 
0.025).  Fig.  5(b)  shows  the  changes  for  Singular  Entropy  Increment 
and  its  differential,  while  the  “o”  refers  to  the  modal  order  identi¬ 
fied.  It  is  easy  to  see  that  Singular  Entropy  Matrix  Pencil  method 
can  get  the  modal  order  exactly. 

When  <7  =  0.025,  we  can  identify  that  the  modal  order  of  signal 
is  6  by  Singular  Entropy  Matrix  Pencil,  and  the  signal  noise  ratio  for 
fitting  is  SNR  =  53.51  dB.  Prony  algorithm  with  actual  modal  order 


is  used  as  a  comparison,  while  its  SNR  is  9.49  dB.  A  comparison  of 
identification  results  is  shown  in  Table  3  and  Fig.  7.  Table  3  shows 
that  the  identified  results  of  Prony  algorithm  are  still  poor  due  to 
its  sensitivity  to  noise  even  if  the  correct  modal  order  is  used. 

From  the  comparison  of  Table  3  and  Fig.  7,  since  Singular 
Entropy  Matrix  Pencil  can  exclude  false  modes,  its  identifica¬ 
tion  results  are  very  close  to  the  actual  signal.  From  Fig.  5,  it 
is  clear  that  the  Singular  Entropy  Increment  of  the  belt  in  case 
2  become  a  little  worse  compared  with  case  1  because  of  the 
noise.  Although,  two  methods  can  extract  the  modal  parameters 
exactly  in  case  1  without  noise.  In  case  2,  the  Singular  Entropy 
Matrix  Pencil  method  performances  well  while  Prony  algorithm 
is  failed. 

3.3.2.  Applied  Singular  Entropy  Matrix  Pencil  method  into 
microgrid 

A  microgrid  simulation  model  is  established  based  on  the  micro¬ 
grid  structure  shown  in  Fig.  1  and  the  proposed  control  strategy  in 


Table  3 

A  comparison  of  identification  results  for  raw  data  with  noise. 


Frequency  (Hz) 

Damping  rate 

Amplitude  (p.u.) 

Phase  (rad) 

Real 

Identified 

Real 

Identified 

Real 

Identified 

Real 

Identified 

Prony  algorithm 

1  20 

20 

0.0036 

0.0046 

1 

1.13 

1.571 

0.181 

2  2 

2 

0.0278 

0.223 

1 

1.91 

1.571 

1.429 

3  0.5 

- 

0.0793 

- 

0.2 

- 

1.571 

- 

Singular  Entropy  Matrix  Pencil 

1  20 

20 

0.0036 

0.0036 

1 

1 

1.571 

1.568 

2  2 

2 

0.0278 

0.0278 

1 

0.99 

1.571 

1.568 

3  0.5 

0.5 

0.0793 

0.0752 

0.2 

0.19 

1.571 

1.564 
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Fig.  8.  The  power  outputs  of  each  DG  and  utility  grid,  (a)  Active  power;  (b)  reactive 
power. 


Fig.  3.  The  parameters  of  the  microgrid  are  shown  in  Table  1 .  In  this 
part,  Singular  Entropy  Matrix  Pencil  method  is  applied  to  analyze 
its  SSS.  Supposed  that  the  microgrid  is  in  interconnected  operated 
mode,  and  the  results  in  isolated  mode  can  be  obtained  analogously. 
The  small  disturbance  is  set  that  the  DGi  has  a  5%  active  power  step 
disturbance  at  10  s. 

Active  power  signals  are  chosen  for  modal  extraction  because 
they  reflect  the  dynamic  of  voltage,  current  and  power  angle 
directly.  The  power  outputs  of  each  DG  are  shown  in  Fig.  8.  We 
extract  the  modal  parameters  of  the  system  by  each  DG’s  active 
power  signal  Pf  (z  =  l,  2,  3)  and  utility  grid  power  Pg  in  the  time 
domain  10-18  s.  The  sample  time  is  Ts  =  0.01  s.  The  identification 


Table  4 

The  main  modes  of  each  DG  extracted  by  Singular  Entropy  Matrix  Pencil. 


Active  power 

SNR (dB) 

Frequency  (Hz)/damping  ratio 

Pi 

168.24 

0.318/0.107,  17.681/0.231 

P2 

155.00 

0.355/0.020,  5.720/0.636,  17.767/0.184 

Ps 

147.457 

0.356/0.016,  5.658/0.621,  16.944/0.209 

Pg 

120.242 

0.327/0.051,  0.594/0.049,  18.261/0.513 

results  are  shown  in  Figs.  9-12.  And  the  main  oscillation  modes  for 
each  signal  are  shown  in  Table  4. 

According  to  Table  4,  it  can  be  observed  that  the  damping  ratio 
of  each  mode  is  positive,  that  is,  the  system  has  no  positive  real 
eigenvalue  and  the  system  is  stable.  It  is  worth  to  notice  that  each 
DG  has  high  frequency  oscillation  modes  whose  frequency  is  larger 
than  5  Hz.  These  modes  reflect  the  rapid  transient  after  system 
disturbance.  On  the  other  hand,  each  DG  has  low  frequency  oscilla¬ 
tion  modes  whose  frequency  is  less  than  5  Hz.  Notation  that  these 
damping  ratios  corresponding  to  low  frequencies  are  small.  Thus, 
these  modes  may  be  more  dangerous  and  easy  to  be  unstable  when 
system  operation  status  changes.  The  DGs  have  different  main  oscil¬ 
lation  modes,  and  that  cannot  be  obtained  by  eigenvalue  analysis. 

4.  Discussion 

The  stability  of  microgrids  is  the  basic  requirement  for  a  safe  and 
reliable  distributed  power  delivery  system.  Unfortunately,  there 
are  a  few  researches  on  it.  Similar  to  the  transmission  network, 
the  stability  for  microgrids  is  divided  into  small  disturbance  stabil¬ 
ity  and  transient  stability  too.  However,  the  operation  conditions  of 
photovoltaic  arrays  and  wind  turbines  depend  on  meteorological 
conditions  highly.  Their  stochastic,  intermittent  and  time-variant 
characteristics  lead  to  the  power  outputs  of  renewable  sources 
change  frequently.  These  can  be  considered  as  small  disturbances. 
Visibly,  small  disturbances  are  so  frequent  that  the  SSS  is  very  com¬ 
plex  and  important  in  microgrids. 

Conventional  SSS  approaches  mainly  use  the  local  linearized 
state  space  model  of  the  nonlinear  system.  Then  the  SSS  results 
can  be  obtained  by  eigenvalues  analysis.  If  all  the  eigenvalues  are 
scattered  beside  left  coordinate  in  the  complex  plane,  the  system  is 
stable.  On  the  contrary,  if  there  has  a  positive  real  part  eigenvalue, 
then  the  system  is  unstable.  However,  due  to  the  unique  charac¬ 
teristics  of  photovoltaic  arrays  and  wind  turbines,  the  linearized 
model  is  time-variant  frequently  even  hard  to  get.  Therefore,  such 
an  approach  is  greatly  flawed  for  SSS  of  microgrids.  In  the  other 
word,  a  further  work  is  strongly  demand  and  a  new  field  need  to 
be  explored  for  SSS  of  microgrids. 

In  this  paper,  a  new  method  for  SSS  of  microgrids  is  pro¬ 
posed,  which  can  extract  the  modal  parameters  from  the  observed 
raw  data.  It  is  verified  by  some  simple  examples  and  a  microgrid 
with  improved  droop  control  strategy.  According  to  the  proposed 
approach,  the  system  is  unstable,  if  there  is  any  mode  with  negative 


x  10'3  x  104 


Fig.  9.  Modal  parameters  extracted  using  active  power  of  DGI  (PI ).  (a)  Singular  Entropy  Increment  and  its  differential  to  get  the  modal  order,  (b)  Comparison  between  raw 
data  and  fitting  data  by  Singular  Entropy  Matrix  Pencil  method,  (c)  Error  for  Singular  Entropy  Matrix  Pencil  method  (R-F).  R,  raw  data;  F,  fitting  data. 
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Fig.  10.  Modal  parameters  extracted  using  active  power  of  DG2  (P2).  (a)  Singular  Entropy  Increment  and  its  differential  to  get  the  modal  order,  (b)  Comparison  between  raw 
data  and  fitting  data  by  Singular  Entropy  Matrix  Pencil  method,  (c)  Error  for  Singular  Entropy  Matrix  Pencil  method  (R-F).  R,  raw  data;  F,  fitting  data. 


Fig.  11.  Modal  parameters  extracted  using  active  power  of  DG3  (P3).  (a)  Singular  Entropy  Increment  and  its  differential  to  get  the  modal  order,  (b)  Comparison  between  raw 
data  and  fitting  data  by  Singular  Entropy  Matrix  Pencil  method,  (c)  Error  for  Singular  Entropy  Matrix  Pencil  method  (R-F).  R,  raw  data;  F,  fitting  data. 


Fig.  12.  Modal  parameters  extracted  using  active  power  of  utility  grid  side  (Pg).  (a)  Singular  Entropy  Increment  and  its  differential  to  get  the  modal  order,  (b)  Comparison 
between  raw  data  and  fitting  data  by  Singular  Entropy  Matrix  Pencil  method,  (c)  Error  for  Singular  Entropy  Matrix  Pencil  method  (R-F).  R,  raw  data;  F,  fitting  data. 


Fig.  13.  An  application  framework  integrated  the  proposed  approach  in  microgrids. 


damping  which  corresponds  a  positive  real  part  eigenvalue.  Thus, 
this  approach  can  be  applied  to  extract  the  SSS  modal  information 
real-time  and  online,  warn  and  emergency  control  in  microgrids. 
The  framework  is  shown  in  Fig.  13  and  the  processes  can  be 
described  as  follows: 

1.  Required  active  power  data  be  measured  in  the  terminal  of  DGs 
firstly.  Then,  the  raw  data  is  transmit  to  microgrid  control  center 
(MGCC). 

2.  After  MGCC  obtains  a  certain  number  of  raw  data  points,  modal 
parameters  can  be  extracted  by  the  proposed  approach. 

3.  According  to  the  modal  parameters,  if  there  is  any  mode  with 
negative  damping  ratio  or  the  damping  ratio  is  less  than  a  pos¬ 
itive  threshold,  the  microgrid  is  judged  as  unstable  or  near  to 
unstable.  Thus,  some  emergency  control  strategies  should  be 
taken  or  warning  should  be  given  according  to  the  contingency 
plans. 

5.  Conclusions 

In  microgrids,  the  SSS  should  attract  more  attention  because 
of  the  unique  characteristics  of  photovoltaic  arrays  and  wind 
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turbines.  This  paper  gives  a  review  on  SSS  of  microgrids  and 
proposes  a  new  SSS  method  based  on  Singular  Entropy  and 
Matrix  Pencil.  Furthermore,  a  microgrid  simulation  model  with  an 
improved  control  strategy  is  used  to  verify  it.  Some  conclusions 
could  be  summarized  as  follows: 

1.  SSS  of  microgrid  is  a  very  important  issue  of  microgrids  due  to 
the  unique  characteristics  of  renewable  sources,  while  there  are 
rarely  researches  on  it. 

2.  Singular  Entropy  Matrix  Pencil  method  uses  the  raw  data 
acquired  from  microgrids.  So  it  can  get  its  real-time  and  online 
SSS  information  dynamically. 

3.  The  proposed  method  applies  Singular  Entropy  theory  to  extract 
the  actual  modal  order  exactly.  It  can  exclude  false  modes, 
reduce  the  computational  complexity  and  improve  the  accuracy. 
Besides,  it  combines  with  Matrix  Pencil  which  can  overcome  the 
signal  noise  very  well  and  extract  modal  parameters  accurately. 

4.  Simulation  results  have  proved  the  feasibility,  validity  and  accu¬ 
racy  of  the  presented  method.  It  shows  that  the  proposed 
approach  is  suitable  for  SSS  of  microgrids  and  the  signal  noise 
ratio  between  fitting  data  and  the  raw  data  is  high. 

5.  A  framework  of  the  proposed  approach  for  modal  parame¬ 
ters  extraction,  monitoring,  warning,  and  emergency  control  in 
microgrids  is  given. 
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